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An analytical method within the frame of linear stability theory is presented for the normal field 
instability in magnetic fluids. It allows to calculate the maximal growth rate and the corresponding 
wave number for any combination of thickness and viscosity of the fluid. Applying this method to 
magnetic fluids of finite depth, these results are quantitatively compared to the wave number of the 
transient pattern observed experimentally after a jump-like increase of the field. The wave number 
grows linearly with increasing induction where the theoretical and the experimental data agree well. 
Thereby a long-standing controversy about the behaviour of the wave number above the critical 
' magnetic field is tackled. 
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PACS numbers: 47.20.Ma, 75.50.Mm 



■ I. INTRODUCTION 

£SJ \ Spontaneous pattern formation from a homogeneous ground state has been studied extensively in many nonlinear 
^ ■ dissipative systems. Among these systems magnetic fluids have experienced a renewed interest in recent years due 
to their technological importance jjj. The most striking phenomenon of pattern formation in magnetic fluids is the 
Rosensweig or normal field instability |^-^| . Above a threshold B c of the induction the initially flat surface exhibits 
a stationary hexagonal pattern of peaks. Typically, patterns are characterized by a wave vector q whose absolute 
value gives the wave number q = |q|. In contrast to many other systems, a comprehensive quantitative theoretical 
and experimental analysis of the dependence of the wave number on the strength of the magnetic field is lacking for 
Q\ \ the normal field instability. There are few but contradictory experimental observations. In experiments where the 
' | • field is increased continuously, there are reports about constant |2||(| as well as about varying wave numbers [[tJ as the 
induction is increased beyond the critical value B c . Notably, all these observations are of entirely qualitative character 

A first theoretical analysis leading to constant wave numbers of maximal growth was presented in |9| . The general 
dispersion relation for surface waves on a magnetic fluid of infinite thickness was analysed for two asymptotic regimes: 
Oh 1 for the inviscid regime and for the viscous-dominated regime. The main result for the latter regime was that taking 
into account viscous effects the wave number of maximal growth is the same at and beyond the critical induction. As 
will be shown below, this argument is rather misleading because realistic fluid properties are not covered by such an 
asymptotic analysis. The two asymptotic regimes in || were combined with very thin as well as very thick layers of 
magnetic fluid and the resulting four regimes were analysed in JTl| . In three regimes a non-constant wave number of 
maximal growth was found. 

All qualitative observations in dUQ] are referring to the final arrangement of peaks. The final stable pattern, 
resulting from nonlinear interactions, does not generally correspond to the most unstable linear pattern. Such a 
pattern should grow with the maximal growth rate and should display the corresponding wave number. Since both 
quantities are calculated by the linear theory, the most unstable linear pattern has to be detected and measured 
experimentally for a meaningful comparison between theory and experiment. No measurements of the most linear 
unstable pattern have yet been undertaken. 

Motivated by this puzzling situation, the paper presents a quantitative theoretical analysis of the wave number 
with maximal growth rate for any combination of fluid parameters. Experimental measurements of the most linear 
unstable pattern are conducted and the data compared with the theoretical results. The system and the relevant 
equations of the problem are displayed in the next Section. Based on the dispersion relation from a linear stability 
analysis, an analytical method is presented to calculate q m and the maximal growth rate w m for any combination 
of material parameters. The details of the method are explained for a magnetic fluid of infinite thickness and the 
results are compared with previous asymptotic results (Sec. III). The method is also applied to magnetic fluids of 
finite thickness (Sec. IV), which allows a quantitative comparison with the experimental data (Sec. V). In the final 
section the results are summarized and further prospects are outlined. 
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II. SYSTEM AND EQUATIONS OF THE PROBLEM 



A horizontally unbounded layer of an incompressible, nonconducting, and viscous magnetic fluid of thickness h and 
constant density p is considered. The fluid is bounded from below (z = —h) by the bottom of a container made of a 
magnetically impermeable material and has a free surface described by z = C,(x,y,t) with air above. The electrically 
insulating fluid justifies the stationary form of the Maxwell equations which reduce to the Laplace equation for the 
magnetic potentials $^ in each of the three different regions. (Upper indices denote the considered medium: 1 air, 2 
magnetic fluid, and 3 container.) It is assumed that the magnetization M^ 2 ) of the magnetic fluid depends linearly on 
the applied magnetic field H^ 2 \ M 1 - 2 ) = (p r — l)H < - 2 - ) , where p r is the relative permeability of the fluid. The system 
is governed by the equation of continuity and the Navier-Stokes equations for the magnetic fluid 

divv = 0, (2.1) 

d t v + (vgrad) v = -- grad (p + p s ) + z/Av + g , (2.2) 
P 

and the Laplace equation in each medium 

A$ w = . (2.3) 

The quantities without an upper index are referring to the magnetic fluid with the velocity field v = (u,v,w), the 
kinematic viscosity is, the pressure p and the acceleration due to gravity g. The first three terms on the right-hand 
side of Eq. ( |2.2| ) result from (1/ p)divT^ where the components of the stress tensor read || 

-p - Mo J (M - pd p M) dH' - p Q — \ S tj + H t Bj +pv{d i v j + djVt) . (2.4) 

The magnetostrictive pressure is given by p s = —po J H pd p MdH' . M, H, and B denote the absolute value of the 
magnetization, the magnetic field and the induction B in the fluid. The governing equations have to be supplemented 
by the appropriate boundary conditions which are the continuity of the normal (tangential) component of the induction 
(magnetic field) at the top and bottom interface 

n (B« -B< 2 )) =0 nxfHW-H^O at z = C (2.5) 

n- (b< 2 ) -B( 3 ') =0 nx (h< 2 ' - H< 3 >) =0 tAz = -h, (2.6) 

the no-slip condition for the velocity at the bottom of the container 

v = d z w = at z = —h , (2-7) 

the kinematic boundary condition at the free surface 

w = <9 t C+ (vgrad) C at z = ( , (2.8) 

and the continuity of the stress tensor across the free surface 

T^) = -aKrij at z = ( . (2.9) 
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The surface tension between the magnetic fluid and air is denoted by a, the curvature of the surface by K = div n, 
and the unit vector normal to the surface by 

n = grad[z-C(z,y,t)] = (-d x t,-d y £,l) 

|grad[*-C(z,y,t)]| V 1 + (dxC) 2 + (dyO 2 ' 



Since the density of air can be neglected with respect to the density of the mag netic fluid and = holds, Eq. ( gj ) 
reduces to 



+ p + po MdH' + p s + y M 1 1 - p v n t {diVj + djVi} = ~oKn 3 , (2.11) 



2 



where is the atmospheric pressure above the fluid layer. In a linear stability analysis all small disturbances from 
the basic state are analysed into normal modes, i.e., they are proportional to exp[— i(ut — qr)]. If 3(o>) > 0, initially 
small undulations will grow exponentially and the originally horizontal surface is unstable. Due to this relation it has 
been established to denote to as growth rate which is in fact true only for its imaginary part. 

Following the standard procedure the linear stability analysis leads to the dispersion relation |Tl|-]T3| (all formulas 
in the references are equivalent to each other) 
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(2.12) 



where fiQ is the permeability of free space, q = y/q 2 — icj /v, and 



_ e^(l + ^)+e-^(l- Mr ) 

_ e^(l + Mr ) 2 - e-*"(l - Mr) 2 ' 1 ' 

The condition of marginal stability, u> = 0, defines the threshold where w chan ges its sign and therefore the normal 
field or Rosensweig instability appears. With u = one obtains from Eq. ( p. 12 ) 

Pg +g 2_ (Mr-l) 2 ^ 2 qA(qh) = Q (2 W) 
(T /io Pr O 

In the limit of an infinitely thick (h — > oo) or an infinitely thin (h — > 0) layer, respectively, the critical inductions are 

r2 _ 2fj, Q fj, r (fj, r + l)y/pag 

Bc >°° ~ ( Pr - l) 2 (2 ' 15) 

wvp^ (2 _ i6) 



whereas in both limits the critical wave number is equal to 

q c = <7c,oo = q C fl = \ — ■ (2-17) 
V a 

The critical values apply to both viscous and inviscid magnetic fluids due to the static character of the instability. 
Based on the dispersion relation ( |2.12| ) the details of the proposed new method are presented exemplarily for a 
magnetic fluid of infinite thickness in the next section. 



III. INFINITE LAYER OF MAGNETIC FLUID 



The starting point of the analysis is the determination of the parameters for which the dispersion relation ( 2.12| ) 
for an infinitely thick layer S 
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(3.1) 



has solutions of purely imaginary growth rates. Such growth rates characterize the viscous-dominated regime described 
by qS 1 0, where 6 — */2vfu) denotes the viscous depth 0. For lu the polar representation u = u)\ + ico2 = 
\uj\ (cos ipo + i sin ipo) is chosen with 
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where n (d) denotes the numerator (denominator) of the argument of arctan. Dimensionless quantities were introduced 
for all lengths, the induction, the time, and the viscosity 



I = q c l 
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where t c is the so-called capillary time. The real and imaginary part of (fO) now read 
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where 



ip = arctan 



\lu\ cos if 



const. 



vq" 



\<jj\ sm ipo 



k = 0,l (3.5) 
fc = 0,l, (3.6) 

(3.7) 



and k distinguishes between the two possible values of the complex root. The value of the constant in ( p.7\ j follows the 
rules of ( |3.2| ). For a purely imagi nary growth rate, uj — iu)2, can take only the two v alue s n/2 (u>2 > 0) and 3n/2 



(uj2 < 0). In the former case Eq. ( p.q ) is always fulfilled, whereas in the latter case Eq. (3.6) holds only if vq 2 > 



For w = iui2 Eq. (3.5) reduces to 



f ± (q, \Q\;D,B) :=[v± 
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(3.8) 



where the ± sign corresponds to <Z>2 ^ 0. The parameters v and B determine the solution of this implicit equation for 
the variables q and \u>\. The solution gives these specific values of the viscosity for which either positive or negative 
purely imaginary growth rates exist (see Fig. [l]). For a supercritical induction of B = 1.05 and q = 1 there exist a 
positive purely imaginary growth rat e fo r all viscosities (fc = 0). Above a critical viscosity, v c = 0.453, a negative 
purely imaginary growth rate solves ( |3.8| ) as well. This critical viscosity also gives the upper bound for the solution 
with k — 1. The value of v c = v c (q, B) increases with increasing induction at constant q and decreases with increasing 
wave vectors at constant B. The critical viscosity is naturally zero at the onset of the instability, v c (q = 1, B — 1) = 0. 

Whereas Fig. [j] shows a situatio n where certain types of solutions of the dispersion relation exist, in Fig. g the 
complete solution of Eqs. (3.5, 3.6) for B = 1.05 and v — 0.037 is plotted (k ~ 0). Around the critical wave number 
q c = 1 a range of wave numbers exist with positive purely imaginary growth rates (filled squares), i.e., there is a band 
of unstable wave vectors. All other growth rates have negative imaginary parts (filled circles, filled diamonds, and 
filled triangles). Therefore the unstable wave vectors are related to the positive purely imaginary growth rates only. 
Focusing on this solution, the imaginary part <I>2 is shown in the vicinity of q c for various strengths of the induction 
and v = 0.037 (v ~ 6.4 • 10 -6 m 2 s _1 ) in Fig. ||. The chosen value of the viscosity characterizes typical magnetic fluids 
in experiments J7|,p^|. All three curves have a maximum in the growth rate u) m = i u)2.m at q-m- One notes that W2,m as 
well as q m are monotonously increasing functions of the strength of the supercritical induction at constant viscosity. 

In order to study the resulting behaviour of 6q = q^/2vj\G)\ at g m , the details of the dependence of ui2, m and q m 
on the field and the viscosity need to be known. The wave number with the maximal growth rate is defined by 
du2/dq — d\oj\/dq = 0. Since \lu\ is given implicitly by f±(q, |cD | ; v, B) — 0, the maximal growth rate results from 
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The cross section of the solutions of Eqs. ( [3.8] , 3.9) gives \u) m \ and q m , which is shown for three different viscosities 
in Fig. Besides a viscosity of real magnetic fluids, two large viscosities (p — 0.4, 2) were chosen to represent the 
regime where the behaviour of the fluid is dominated by the viscosity. For all three viscosities the wave number q m is 
not constant, i.e., for finite viscosities q m depends on the external control parameter B. With increasing viscosity q m 
varies less with increasing induction, e.g. Aq rn — 1.68 for v — 0.037 reduces to A.q m — 0.29 for v — 2 at an induction 
difference of AB = 0.5. For small viscosities q m depends linearly on B if B is not too large. This linear dependence 
shifts towards higher values of B with increasing viscosity (compare v — 0.037 and 0.4). At the largest viscosity, 
v ~ 2, no linear behaviour can be observed for 1 < B < 1.5. 

The analysis reveals that only in the case of infinitely large viscosities (with respect to the viscosity of real magnetic 
fluids) a constant wave vector of maximal growth q rn = 1 can be expected. Taking into account viscous effects does not 
necessarily lead to a constant q m . For a better comparison with [^j the value of qS at q m is calculated and is plotted 
for the three viscosities chosen in Fig. |j| The graphs show clearly that q m 6 ^ 1 holds only in the close vicinity of the 
critical induction for large viscosities (v = 0.4, 2) and in the limit B — 1 (|o) m | = 0, 5 = oo) for realistic viscosities 
(v = 0.037). Because 0.43 < q m S < 1.08 for v = 0.037, realistic fluid properties are neither covered by the other 
asymptotic regime q m S <C 1 analysed in Therefore the experimental observation in [^|^] can not be explained 
by the result of an asymptotic analysis which does not meet the features of the experimental fluids. By plotting 

the know analytical result in the inviscid regime, q m = (1/3) (2B 2 + V4B 4 - 3) pH , realistic magnetic fluids tend 

rather to the limit v — > than to the limit v — > 00 (see Fig. ^) as exploited in earlier studies gj. But for quantitative 
comparisons in typical experimental setups, asymptotic analyses p,p!T[ are insufficient. 



IV. FINITE LAYER OF MAGNETIC FLUID 



Since the experiments are performed with a vessel of finite depth, the method presented in the previous section has 
to be applied to magnetic fluids of finite thickness. With the polar representation of u), the real and imaginary part 
of the dispersion relation (2.12) read 
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The explicit form of the abbreviations Rj^Ii, N±, N 2 , and (i = 1, . . . , 4) is deferred to Appendix A. For purely 
imaginary growth rates, oj = iuo 2l Eq. (4.2) is fulfilled without any restrictions for to 2 > as well as for lo 2 < 
in contrast to the case of an infinitely thick layer. Only positive purely imaginary growth rates are of interest for 
comparison with the experiment. Therefore the function / is now of the form 
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where this implicit equation for q and lo contains the additional parameter h. Fig. ^ shows the solution for three 
different depths of the layer at a supercritical induction of B = 106 • 10~ 4 T. The used material parameters of the 
commercially available magnetic fluid EMG 901 (Ferrofluidics Corporation) are listed in Table |. The graphs show 
that the wave number of maximal growth q m clearly varies less with h than the maximal growth rate uo 2 . m itself. 
Notably, the solution for a layer of 2 mm thickness is already near the infinite case illustrated by h — 100 mm. 
Therefore a filling with h > 1 cm of magnetic fluid can be considered as an infinite thick layer. To make such an 
estimate is an asset in the use of the complete equations. Because they cover the entire range of thickness, < h < 00, 
in extension to the asymptotic analysis h ~ and h ~ 00 in ]ll[ |. 

To analyse the behaviour of u> 2 ^ m and q m on B and h, the maximal growth rate, given by df+/dq = 0, has to be 
determined. As the resulting i mplic it function is quite lengthy thus we do not give the explicit form here. The cross 
section of the solutions of Eq. ( |4.3|) and df + /dq = leads to w 2 , m and q m . Their dependence on the supercritical 
induction and the thickness of the layer is shown in Fig. |?]. The wave number of maximal growth increases linearly 
with B with the exception of B near the height dependent critical value B c ^. The linear behaviour is independent of 
the thickness of the layer and holds up to 30% above B c ^- The maximal growth rate uj 2%m starts to grow like a square 
root above the onset of the instability. This square root behaviour becomes less pronounced with thinner layers. 
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Through the implicit character of the functions an analytical expression cannot be given for the dependence of q m 
and u>2, m on B and h. Alternatively, a two parameter fit is tested, which describes the generic behaviour over a wide 
range of layer thicknesses. An excellent agreement is achieved for h > 4 mm by 

q m = 3.265 - 0.09VB for 0.001 < B < 0.2 (4.4) 
w 2 ,m = 1.18^5 + 2.9B for 0.001 < B < 0.2 , (4.5) 

(see Fig. ||) where B = (B — B c ,h) / Bc,h, Qm = {q m — q c ,h) / q c ,h, and u>2,m = ^2, m t c denote the scaled distances from the 
critical values. For small B the behaviour of q m is only weakly nonlinear whereas the behaviour of u)2,m is determined 
by the square root term. A careful inspection of the data reveals that for h = 2 mm (filled circles) small deviations 
from the proposed fits appear: q m grows linearly over the entire B region (see insert in Fig. (a))- Thus h — 2 mm 



indicates the lower limit of the validity of (4.4, 4.5) 



Since the fit covers the region of infinite thick layers one can expand (3.5, 3.£) for small B, q m , and u>2.m- Taking 
into account that the dimensionless viscosity is also small for real magnetic fluids {p ~ 0.0483 for EMG 901) the 
expansion leads to 

6 18 u 2 

q m = —B 1 -—B i for0<B«— (4.6) 
v 6 



2 - 3 v 2 
-B - -=sB for < B < - 
v v 6 



w 2 , m = —B — rrrB for < B < — , (4.7) 



where the coefficients depend on the viscosity. Fig. ^ shows the very good agreement between the numerical solution 
and the expansion for B <C v 2 /6. The region where the expansion holds extends with the square of the viscosity. The 
expansion (4l| 4.7) and the scaling ( |4.4| , i.l) show that the behavior of q m and u>2, m is entirely governed by the two 



parameters viscosity and induction for not too thin layers. The third parameter h has only little effect in this regime. 

From the results shown in Fig. ^ one notes the height dependence of B c and q c at the onset of the instability (see 
also Fig. 4, 5 in This dependence for B c can be exploited to measure the permeability of the magnetic fluid 

'just in time' for the experiment. Since the quotient of the two limits ( 2.15 , |2.16 ) depends on fi r only 



Bc,0 _ 2^r (4 g) 



Be, 00 l^r 1 

the determination of the two limits of B c offers a feasible access to \i r of the magnetic fluid. From Fig. |Io| it can be 
seen that B c increases monotonously from B CjOC} towards B c q with decreasing layer thickness. Since the preparation 
of very thin layers is laborious and delicate it would be desirable to shift the thin layer limit towards thicker films. 
This can be achieved by an increase of the surface tension. A modified surface tension is accompanied by changes in 
the density and permeability of the fluid. But these changes are of a much smaller scale than those of the surface 
tension. The modified viscosity does not affect the determination of B c h- Therefore the surface tension is changed 
whereas all other quantities remain constant. By increasing the surface tension by a factor of 10 (100), B c ^ may be 
measured for layers nearly 1 (2) orders of magnitude thicker than for a system with the original surface tension (see 
Fig.0). 

Abou et al. analysed the limit of thin films of magnetic fluid for vanishing and infinitely large viscosities. In both 
cases the analytical result Q 

1 / lit i 11 / I,. 1 \2 \ 

(4.9) 

is the same. Therefore, one can assume that the dependence of q m on B is not influenced by the viscosity in the 
thin film limit. Since the present method allows to calculate q m for any combination of parameters, we are able 
to accomplish a test of this assumption. As Fig. [ll] shows, the behaviour of q m on B is indeed independent of v 
(The tested viscosities cover a range from 6.5 • 10~^n 2 s _1 to 2 • 10~ 5 m 2 s -1 .). The variation of q m on the applied 
magnetic induction in thin films was measured in earlier experiments where the magnetic fluid either was prepared at 
the bottom of a quartz chamber [jl6| or was laid on top of a denser fluid Jlj],[l8| ■ In all three experiments the spacing 
between the final arrangement of peaks was measured in dependence of the applied field and a nonlinear behaviour 
was found. 




6 



V. MEASUREMENT, RESULTS, AND COMPARISON WITH THEORY 



In this section we report on experimental results of the dependence of the maximal wave number on the supercritical 
magnetic induction. First we present the experimental setup, next we give a characteristic example of the pattern 
evolution. We continue with a description of the techniques applied to extract the wave number of the patterns. 
Finally the experimental results are compared with the theoretical results of the previous section, particularly the 
predicted growth of the maximal wave number. 

Our experimental setup is shown in Fig. |l2|. A cylindrical Teflon® vessel with a diameter of d = 12 cm and a 
depth of 3 mm is completely filled with magnetic fluid and situated in the center of a pair of Hclmholtz coils. The 
experiments were performed with EMG 909. The fluid is illuminated by 90 red LEDs mounted on a ring of 30 cm 
diameter placed in a distance of 105 cm above the surface. A CCD-camera is positioned at the center of the ring. By 
this construction a flat fluid surface reflects no light into the camera lens, however, an inclined surface of proper angle 
will reflect light into the camera ]l9]j . The CCD-camera is connected via a framegrabber to a Pentium 90 MHz PC 
and serves additionally as a fundamental clock for timing the experiment. In the theoretical analysis the supercritical 
magnetic field is assumed to be instantly present, thus in the experiment the magnetic field has to be increased 
jump-like from a subcritical value B to the desired value B. For all measurements B was fixed to 133 • 10~ 4 T. 
The jump- like increase of the field is initiated by the computer. Its D/A converter is connected via an amplifier 
(fug Elektronik GmbH) to the Hclmholtz coils (Oswald Magnetfeldtechnik) . The magnetic system can not follow the 
control signal instantly, its relaxation time tb to a jump-like increase of the control signal depends on the jump hight 
AB = B — Bq. For a maximal jump hight of AB — 70 • 10~ 4 T the relaxation time mounts up to tb = 80 ms. The 
other characteristic time scales of the system are the capillary time scale, t c ~ 13 ms, and the viscous time scale, 
t v = l/(ql v) ~ 450 ms. 

For the empty Helmholtz-coils the spatial homogeneity of the magnetic field is better than ±1%. This grade is 
valid within a cylinder of 10 cm in diameter and 14 cm in hight oriented symmetrically around the center of the 
coils. Two Hall-probes are positioned immediately under the Teflon® dish. A Siemens Hall-probe (KSY 13) serves to 
measure the magnetic field during its jump-like increase, and is connected to the digital voltmeter (Prema 6001). For 
measuring a constant magnetic field and for calibration purpose we use a commercial Hall-probe (Group3-LPT-231) 
connected to the digital teslameter (DTM 141). Both devices are controlled via IEEE-bus by the computer. 

Next we give a characteristic example for the evolution of the surface pattern during a jump-like increase of the 
magnetic field. Fig. |l3|(a) shows circular surface deformations taken At = 180 ms after the start of the experiment. 
This surface deformations are first created at the edge of the dish, because of the discontinuity of the magnetic 
induction induced by the finite size of the container. The circular deformation is fixed in space, and its amplitude 
grows during the jump-like increase of the magnetic field. With increasing time more circular deformations evolve, 
approaching the center of the dish (see Fig. |l3|(a)). Onto this pattern, Rosensweig peaks emerge out of the crests of 
the circular surface deformation, as can be seen in Fig. |l3|(b). After this transient concentric arrangement a hexagonal 
pattern of Rosensweig peaks evolves (see Fig. |l3|(c)). 

The theoretical results stem from a linear theory which can determine correctly the critical values of the pattern 
selected by the instability only at the threshold. Above the threshold, a band of wave numbers will become unstable, 
where the mode with the largest growth rate is the most unstable linear mode of the flat interface. Due to nonlinear 
effects the final stable pattern docs not generally correspond to the most unstable linear mode, as shown here and 
in other experiments p0| . Therefore it has to be stressed that for the comparison with the linear theory not the 
stable hexagonal pattern, but the most early stage of the pattern, namely the transient circular deformations are 
most appropriate. 

The wave number of the circular deformations is extracted from the pictures in the following way. First the algorithm 
scans the diagonals of the picture for the local maximum of the grey levels. Starting at the corners of the pictures, 
it detects points which are situated at the edge of the Teflon® dish. Two of the edge points are marked by white 
circles in the left part of Fig. |l4|(a). From the full set of four points the algorithm calculates the center of the dish 
denoted by the half circle at the right part of Fig. |l4|(a). In order to control the precision of the algorithm a white 
circle with the proper radius of the dish is constructed around the detected center. In the next step we calculate the 
radial distribution of the grey values for all pixels within the circumference of the dish, as shown in Fig. |lj(b). For 
comparison Fig. 14|(c) gives an artificial two dimensional representation of the grey value distribution. 

In all of the three pictures of Fig. [l4| one can easily discriminate three zones. In the innermost zone only small 
surface undulations exist, which give rise to the unstructured part of the grey value distribution in Fig. [l4](b,c). For 
larger radii one finds the area of circular surface deformations which generates the bi-periodic peak pattern in the 
distribution. The large peaks are correlated with the deformation troughs and the small peaks with the deformation 
crests. Finally the outermost zone includes the edge of the dish together with the first, edge induced, deformation 
crest. For estimation of the wave number we discard the inner- and outermost zones. The top part of the peaks in 
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the remaining zone is fitted by a polynomial of second order. The average distance of their maxima gives the half 
of the desired wavelength. For the rather small number of deformations available, the above presented method of 
wave number calculation turned out to be more stable and to give more precise results than the competing method 
of two dimensional Fourier transformation. Together with the picture the momentarily magnetic induction has been 
recorded during the jump-like increase. This allows an exact relation of the extracted wave number to the instantly 
prevailing magnetic induction. 

Let us now focus on the experimental results displayed in Fig. [l5|, where the wave number q is plotted versus the 
magnetic induction B. Each open square denotes the wave number extracted from a picture taken during a jump-like 
increase of the magnetic field to B > B c . The estimated maximal errors for q of ±4.2% and for B of ±0.9% are not 
plotted for the purpose of clarity. The dashed line displays the theoretical results for the listed material parameters 
of the magnetic fluid EMG 909. Using fj, r as a fit-parameter gives the solid line with \x r ~ 1.86. The fitted value for 
\x r differs by 3.3% from the value given by Ferrofluidics, a deviation which is well within the tolerance of production 
specified by Ferrofluidics. Obviously there is a rather good agreement between the experimental results and the 
theoretical graph for \x r ~ 1.86 marked by the solid line. The linear increase in the appearing wave number, both in 
experiment and in theory, is our main outcome. 

Comparing the two theoretical curves in Fig. an increase of /i r results in a decrease of the critical magnetic 
inducti on, w hereas the critical wave number remains constant, as can bee seen fromEqs. fl2.15j |2.16| , |2~17| ). According 



to Eq. (2.17) a constant critical wave number implies that the density and the surface tension are constant. Therefore 
we refrain to consider them as additional fit parameters. As can be seen from Fig. |] changes in the viscosity by an 
order of magnitude are necessary to cause a relevant influence on the behaviour of wave number of maximal growth 
on the induction. Therefore changes in the viscosity due to small thermal fluctuations in the experiment can be 
neglected. 

We find a linear wave number dependence of the circular surface deformations. This pattern is a more simple 
realization of the normal field instability than the familiar hexagonal pattern of Rosensweig cusps. The latter one 
is obtained by a symmetrical superposition of three patterns of parallel stripes with the wave vectors separated by 
120 degree j21| . Obviously the circular surface deformations can be regarded as a stripe pattern favoured by the 
symmetry of the dish. As a consequence they appear first, before nonlinear interactions select in a later stage the 
hexagonal pattern. This situation is well known from Rayleigh-Benard convection in cylindrical containers, where, due 
to side-wall induced convection, concentric target patterns appear instead of hexagonal structures. Our observations 
agree in part with recent findings by Browaeys et al. ]22^j . They detected circular surface deformations for a constant, 
subcritical magnetic field of 0.797? c . In contrast to their experiment, we do not perform a periodic modulation, but a 
jump-like increase of the magnetic induction. Thus we have no interference with additional waves propagating onto 
the circular deformations. Therefore a measurement of the wave number, as described above, could be realized. 

The circular surface deformations have to be distinguished from circular, meniscus induced surface waves emitted 
from the edge of lateral cell walls p3[ . Here, the circular deformations are induced by the discontinuity of the magnetic 
induction at the edge of the container. The formation of a meniscus is eluded by a brimful filling of the dish and 
by the design of the vessel which has a slope with respect to the horizontal of 15°, the contact angle between the 
magnetic fluid and Teflon®. 

Finite size effects due to the finite size of the vessel are rather small in the experiment. Applying the arguments of 
Edwards and Fauve |24j] , the width of the band of unstable wave numbers, Aq ~ 5.1 cm - 1 for B = 180 • 10" 4 T and 
h = 3 mm, is much larger than ir/d ~ 0.3 cm , the wave number separation between the quantized modes of the 
vessel. Thus the influence of the vessel size can be neglected and the developing pattern is insensitive to the vessel 
size. 

For the experiments we have chosen a magnetic fluid with a rather low value of magnetic permeability fi r = 1.86, 
in order to keep hysteresis effects small. Indeed, with our resolution a hysteresis can not be detected. Hysteresis 
strength proves to increase monotonically with the permeability of the magnetic fluid ]25) . Thus, the influence of 
higher permeability on the wave number of maximal growth remains to be investigated experimentally in the context 
of a more complex situation of a transition with large hysteresis. 



VI. SUMMARY 



By means of the polar representation of the complex frequency lj the dispersion relation for surface waves on viscous 
magnetic fluids is split into a real and an imaginary part. The parameters are determined for which pure imaginary 
solutions lo = iu)2 and ll>2 > for both parts exist. For these parameters the originally horizontal surface is unstable, 
because initially small undulations of the surface, proportional to exp(— iut), grow exponentially. The imaginary part 
of the dispersion relation is fulfilled mostly automatically. From the real part, the wave number with maximal growth 
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rate q m and the maximal growth rate oj 2 ^ m itself can be easily determined. It can be done for any combination of 
material parameters and for any thickness of the layer. This is the strength of the presented analytical method which 
covers the entire parameter space between the previously studied asymptotic cases . It therefore allows to study 
the transition from one limit to the other. Such a transition is exemplarily illustrated for an infinitely thick layer with 
viscosities varying between zero and infinity. 

For magnetic fluids of infinite depth it is shown that earlier qualitative observations of constant wave numbers 
above the critical magnetic field [||](| cannot be explained by the result of an asymptotic analysis Q . The analysis 
in H does not cover the features of the experimental fluids. In order to apply a theory, where the field is instantly 
present, a jump-like increase of the field in the experiments is essential. Therefore the results for a continuously 
increased field m are inappropriate for a comparison with such a theory. Furthermore, we have demonstrated that 
the transient pattern is the most suitable one to be compared to the linear theory. Taking all these into account, we 
are able to observe a linear increase of the wave number of maximal growth with increasing magnetic induction. This 
linear increase is quantitatively confirmed by the linear theory. 

An increasing wave number with increasing field was also observed in the corresponding electrical setup where a 
liquid metal is subjected to a normal electric field but this result is based only on a qualitative observation. The 
authors emphasize as well the importance of a fast build up of the field. 

It is very attractive to test in further experiments whether the predicted generic behaviour of the maximal growth 
rate can be confirmed in the weakly nonlinear regime. As an expected outcome, u> m ,exp should start to grow like a 
square root with increasing supercritical induction. Furthermore, it remains to be seen whether the linear increase 
of the wave number of the linearly most unstable pattern lasts in the final hexagonal pattern. A confirmation would 
mean that also the wave number of the final pattern varies if the induction is jump-like increased. We point out, 
that for a continuous increase of the induction the behaviour of the wave number for both the transient and the final 
pattern remains to be elucidated. 
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VII. APPENDIX 



The abbreviations in Eqs. (4.1 4.2) read explicitly 



Ri 
h 
Ni 
R 2 



N 2 

h 

h 



[cosh(2<?i/i) — cos(2q 2 ^)] {<?i sinh(2gi/i) + q 2 sin(2q 2 h) — qcoth(qh) [cosh(2qi/i) — cos(2q 2 h)]} , 
[cosh(2<7i/i) — cos(2<72^)] [<?2 smh(2q\h) — q~i sin(2q2^)] , 

{qi sinh(2<7i/i) + q 2 s\xi{2q 2 }i) — qcoth(gft.) [cosh(2qi/i) — cos(2q 2 h)]} 2 + {q 2 sinh(2qi/i) — q\ sin(2q2^)} 2 



6g 2 <?2 



[q~l sinh(2gi/i) + q 2 sin(2q 2 h)] 



Q1Q2 ~ 4<zi<22] sinh(2<7i/i) - qi su\(2q 2 h)] 



[5q 2 + 2q 2 {q\ 
- [Aq 2 qxq 2 + 

[iq 2 q!q 2 + iqfq 2 - 4giqf] [q x sinh(2qi/i) + q 2 sin(2q 2 h)] 
+ [5q 2 + 2q 2 {q 2 - q 2 ) + qf - 6q 2 q 2 + <?|] [q 2 sinh(2g 1 / l ) 
cosh(2<7i h) — cos(2q 2 h) , 
qtanh{qh) [6q 2 (q 2 - q 2 ) + q A + q\ - 6q 2 q 2 
qt&nh{qh) [Uq 2 ^ + ±q\q 2 - 4giqf] , 
Aq 2 sinh(gift) cos(q 2 h) [qi(q 2 + q\ - q%) 
Aq 2 sinh(gi/i) cos(q 2 h) [q 2 (q 2 + q\ - 
cosh(qh) [cosh 2 (q\h) — cos 2 {q 2 h)\ . 

In ( |7TT| - f7. 1 1| ) the shorthands 



<Z2 4 ] 



e 2 ) 



2q x q 2 ] 
2q 2 g 2 ] 



qi sm(2q 2 h)} 



Aq 2 cosh(qih) sm(q 2 h) [q 2 (q 2 
Aq 2 cc-sh^h) s'm(q 2 h) [qi(q 2 



qf) + 2q 2 q 2 ] 

ql) - 2?igl] 



(7.1) 
(7.2) 
(7.3) 

(7.4) 

(7.5) 
(7.6) 
(7.7) 
(7.8) 
(7.9) 
(7.10) 
(7.11) 



9i = 




uj\ 2 cos 2 ipo fib + 2kir 
— cos 1 



(7.12) 







4 f 2 |w|sin^ \ . |w| 2 cos 2 </> . / V + 2/ot x . 
q-2 = \ [q 2 + + -o sm (7.13) 



were used where 



, — \0J COSCA) /_ ,n 

ip = arctan — : — — h const . ( 7 .14) 

q z v + \lu\ sm ip 
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FIG. 1. Purely imaginary growth rate Co = iu)2 as function 
of the viscosity for q = 1. For k — (solid lines) positive 
values of u)2 exist for all viscosities whereas negative ones 
exist only above a critical viscosity v c = 0.453. For k = 1 
(dashed line) the upper bound for Co — iu)2 is given by v c . 
The long-dashed line indicates the condition uq 2 = for 
u> 2 < 0. 




0.0 0.5 1.0 1.5 2.0 



wave number q (c^) 

FIG. 2. Growth rate uj — u>\ + i u)2 as function of the wave 
number q for v — 0.037 and B = 1.05. The empty symbols 
(a, o, o) show 5R(w) = u>i and the filled symbols (A, ♦, •, ■) 
display Q(co) = UJ2. A positive purely imaginary growth rate 
(■) exists only in the vicinity of q = 1. 
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1.0 1.5 2.0 

wave number q (q c ) 

FIG. 3. Positive purely imaginary growth rates as a func- 
tion of the wave number q and the strength of the induction 
B for v — 0.037, typically for magnetic fluids in experiments 
[^,^|. The maximum is given by ui2,m and q m which both 
increase monotonously with B. 
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FIG. 4. Maximal wave number q m as function of the su- 
percritical induction B for different viscosities. q m is a 
monotonously increasing function of B with the exception 
Qm — 1 §1 in the case of infinitely large viscosities (lower 
dot-dashed line). In the limit of an inviscid fluid (upper 
dot-dashed line) the dependence of q m on B is given by 

q m = (1/3) (2B 2 + \/4B 4 - 3 ) ju). 
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induction B (B ) 

v c,oo' 

FIG. 5. Behaviour of qd at q m as a function of the induc- 
tion for the viscosities v = 0.037, 0.4, and 2. The condition 
QmS ^> 1 holds only in the close vicinity of the critical induc- 
tion, B w 1, for large viscosities (p = 0.4, 2) and in the limit 
B — 1 for realistic viscosities (u = 0.037). 
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FIG. 6. Positive purely imaginary growth rates as function 
of the wave number q and the thickness of the layer at a 
constant induction of B = 106- 10 -4 T. Remarkably the graph 
for a layer of 2 mm is already close to the limit of an infinite 
thick layer illustrated by h = 100 mm. Whereas the wave 
number of maximal growth shows only a small variation, the 
maximal growth rate itself displays more distinct changes. 
Material parameters of the fluid EMG 901 are listed in Table 
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FIG. 7. Maximal wave number q m (a) and maximal growth 
rate u>2,m (b) as function of the supercritical induction 
B > B Cy h for three different thicknesses, (a) q m increases 
linearly with B except for B near the critical value B Ct h- The 
area of nonlinear behaviour shrinks with the shrinking thick- 
ness of the layer, (b) U2,m starts to grow like a square root 
above the onset of the instability. This square root behaviour 
becomes less pronounced with thinner layers. Material pa- 
rameters of the fluid EMG 901 are listed in Table || 
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scaled induction B 

FIG. 8. Scaled maximal wave number g m (a) and scaled 
maximal growth rate d>2,m (b) as function of the scaled su- 
percritical induction B. The data are calculated for h — 100 
mm (o), 50 mm (*), 10 mm (+), 4 mm (□), 2 mm (•). For 

h > 4 mm the data are fitted by q m = 3.26B - OmV B for 

q m (solid line (a)) and by &2,m = 1.18v5 + 2.9B for il>2,m 
(solid line (b)). Small deviations from the generic behaviour 
can be seen for h — 2 mm (insert). Material parameters of 
the fluid EMG 901 are listed in Table @. 
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FIG. 9. Scaled maximal wave number g m (a) and scaled 
maximal growth rate o>2,m (b) as function of the scaled su- 
percritical induction B. The data are calculated for h — 100 
mm (o), 50 mm (*), and 10 mm (+). For B < v' 1 /6, i.e. 
B < 4 • 10~ 4 for EMG 901, the agreement with the analyt- 



ical results q m = (6/v 2 )B 2 - (18/^ 4 )B' i (solid line (a)) and 
^2,m = (2/v)B - {3/v 3 )B 2 (solid line (b)) is very good. Ma- 
terial parameters of the fluid EMG 901 are listed in Table 
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FIG. 10. The critical induction B c> h versus the thick- 
ness of the layer h for three different surface tensions: 
a = 2.275 • 10~ 2 kgs~ 2 (solid line), a = 2.275 • lO^kgs" 2 
(dashed line), and a = 2.275 kg s -2 (long-dashed line). By 
increasing the surface tension by a factor of 10 (100), -B c ,o 
can be measured for layers nearly 1 (2) orders of magnitude 
thicker than for a system with the original surface tension. 
The remaining material parameters of the fluid EMG 901 are 
listed in Table EL 
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FIG. 11. Maximal wave number q m as function of the 
supercritical induction B of a thin film, h — lfim, for 
three different viscosities: v = 6.5 • 10 _8 m 2 s _1 (o), 
v = 6.5 • 10~ 6 m 2 s _1 (□), and v = 2.0 ■ 10 -5 m 2 s^ 1 (A). 
The numerical data show that the behaviour of q m on B is 
independent of v and it is given by <jm = (1/4) (c+V^S) 
with c = 3B 2 (n,. + l)/(2/Lt r ) |ll[ (solid line). The remaining 
material parameters of the fluid EMG 901 are listed in Table 
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FIG. 12. Scheme of the experimental setup 




FIG. 13. Series of snap-shots of the principal pattern evolution of the magnetic fluid for a jump from Bo < B c to B > B c , 
illuminated from above by a ring of LEDs. The pictures are taken At = 180 ms (a), 280 ms (b), 560 ms (c) after the start of 
the increase of the magnetic field. 




FIG. 14. Three steps of the picture processing to extract 
the wave number: (a) Reflections of circular surface defor- 
mations. The white line in the left upper and lower corner 
marks the calculated edge of the dish, the circles on the edge 
serve to calculate the center of the dish, (b) Radial grey level 
distribution of (a), (c) Two dimensional representation of (b). 
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FIG. 15. Plot of the wave number q versus the magnetic 
induction B. The open squares give the experimental values, 
the dashed line displays the theoretical results for the material 
parameters of EMG 909 listed in Table [[]. Using fi r as a 
fit-parameter gives the solid line. 



TABLE I. Material parameters of EMG 901 and EMG 909. 



EMG 901 Source EMG 909 Source 

fi r 4.0 Ferrofluidics 1.8 Ferrofluidics 

p [kg-m -3 ] 1.53-10 3 [[| 1.02-10 3 Ferrofluidics 

v [m 2 ^ 1 ] 6.54-10" 6 Ferrofluidics 5.88-10" 6 Ferrofluidics 

a [kg-s~ 2 ] 2.27-10" 2 || 2.65-10~ 2 @ 
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